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ABSTRACT Results are presented supporting a regulatory role for the product of the MA3302 gene locus (designated MreA) previ- 
ously annotated as a hypothetical protein in the methanogenic species Methanosarcina acetivorans of the domain Archaea. Se- 
quence analysis of MreA revealed identity to the TrmB family of transcription factors, albeit the sequence is lacking the sensor 
domain analogous to TrmBL2, abundant in nonmethanogenic species of the domain Archaea. Transcription of mreA was highly 
upregulated during growth on acetate versus methylotrophic substrates, and an mreA deletion (AmreA) strain was impaired for 
growth with acetate in contrast to normal growth with methylotrophic substrates. Transcriptional profiling of acetate-grown 
cells identified 280 genes with altered expression in the AmreA strain versus the wild-type strain. Expression of genes unique to 
the acetate pathway decreased whereas expression of genes unique to methylotrophic metabolism increased in the AmreA strain 
relative to the wild type, results indicative of a dual role for MreA in either the direct or indirect activation of acetate- specific 
genes and repression of methylotrophic-specific genes. Gel shift experiments revealed specific binding of MreA to promoter re- 
gions of regulated genes. Homologs of MreA were identified in M. acetivorans and other Methanosarcina species for which ex- 
pression patterns indicate roles in regulating methylotrophic pathways. 

IMPORTANCE Species in the domain Archaea utilize basal transcription machinery resembling that of the domain Eukarya, rais- 
ing questions addressing the role of numerous putative transcription factors identified in sequenced archaeal genomes. Species 
in the genus Methanosarcina are ideally suited for investigating principles of archaeal transcription through analysis of the ca- 
pacity to utilize a diversity of substrates for growth and methanogenesis. Methanosarcina species switch pathways in response to 
the most energetically favorable substrate, metabolizing methylotrophic substrates in preference to acetate marked by substan- 
tial regulation of gene expression. Although conversion of the methyl group of acetate accounts for most of the methane pro- 
duced in Earth's biosphere, no proteins involved in the regulation of genes in the acetate pathway have been reported. The re- 
sults presented here establish that MreA participates in the global regulation of diverse methanogenic pathways in the genus 
Methanosarcina. Finally, the results contribute to a broader understanding of transcriptional regulation in the domain Archaea. 
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Microbes in the domain Archaea utilize basal transcription 
machinery resembling that of members of the domain Eu- 
karya which includes a TATA box promoter sequence, a TATA 
box-binding protein (TBP), a homologue of the transcription 
factor TFIIB (TFB), and a polymerase (Pol) Il-like RNA polymer- 
ase containing between 8 and 13 subunits (1). In contrast, the 
majority of characterized and putative transcription factors in the 
domain Archaea are homologous to activators and repressors in 
the domain Bacteria (1-8). A recent analysis of 52 sequenced ar- 
chaeal genomes revealed 3,918 putative TFs, among which the 
HTH_3, AsnC, TrmB, and ArsR families were universally abun- 
dant in all 52 genomes (9). Interestingly, the proportion of small 
TF's (100 to 200 residues) was considerably larger than the pro- 
portion encoded in bacterial genomes. This understanding raises 
basic questions regarding the physiological roles of putative TF's 
and their transcriptional regulatory networks. Although 



methane-producing species (methanogens) are the largest and 
most broadly studied group in the domain Archaea, relatively few 
are archaean specific (8, 10, 11, 12), and bacterium-like (13, 14) 
TF's have been identified which regulate genes involved in elec- 
tron transport (13, 14), nitrogen fixation and assimilation (8, 11), 
oxidative stress (10), and biosynthesis of tryptophan (12). 

Methanogens have ancient origins and are terminal organisms 
of anaerobic microbial food chains which function in diverse an- 
aerobic environments of Earth's biosphere to convert complex 
organic matter to methane, an integral component of the global 
carbon cycle. The nearly 1 billion metric tons produced annually 
contributes to the accumulation of atmospheric methane, a green- 
house gas severalfold more effective than C0 2 in radiating light 
energy back to Earth (15). Implementation of this natural process 
for the large-scale conversion of renewable biomass and organic 
wastes to methane is a viable alternative to the use of fossil fuels. 
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Microbes, primarily of the domain Bacteria, initiate the food chain 
by decomposing complex organic matter to acetate, methyl- 
amines, methanol, formate, H 2 , and CO z , which are growth sub- 
strates for methanogens. Species of the genus Methanosarcina are 
the only methanogens known to grow and produce methane from 
acetate and the methylotrophic substrates methanol, methyl- 
amines, and methylsulfide. These versatile methanogens switch 
pathways in response to availability of the most energetically fa- 
vorable substrate, metabolizing methylotrophic substrates in 
preference to acetate, as marked by substantial regulation of gene 
expression (16, 17). Thus, Methanosarcina species are ideally 
suited for investigating complex regulation in the domain Ar- 
chaea. Furthermore, robust genetic systems are available for sev- 
eral Methanosarcina species (18-20). 

Although much is known of enzymes essential to the diverse 
methanogenic pathways of Methanosarcina species, relatively little 
is known regarding regulation of gene expression. A recent genetic 
investigation identified a class of regulatory proteins in Methano- 
sarcina acetivorans which coordinate expression of methanol- 
grown and methylsulfide-specific methyltransferase genes (21, 
22). However, although acetate is the major substrate for biolog- 
ical methane formation in Earth's biosphere, no proteins involved 
in the regulation of genes essential to the pathway for methano- 
genesis from acetate have been reported. 

A previous proteomics investigation of M. acetivorans C2A 
identified numerous proteins with elevated levels during growth 
on acetate versus methanol, indicating potential roles during 
growth on acetate (17). Bioinformatics analysis of one of the pro- 
teins (encoded by the MA3302 gene locus), previously annotated 
as hypothetical and upregulated 250-fold, revealed sequence iden- 
tity to the TrmB family of archaeal regulatory proteins. Thus, 
MA3302 was investigated with a combination of genetics, com- 
parative bioinformatics, and transcriptomics to ascertain the po- 
tential for regulation of genes unique to acetate metabolism and to 
broaden the understanding of regulation in the domain Archaea. 
The results reveal that the product of MA3302 is a previously 
unrecognized transcriptional regulatory protein that we designate 
MreA (Methanosarcina regulator of energy-converting metabo- 
lism). MreA was shown to be involved in the activation of genes 
encoding enzymes unique to the pathway of methanogenesis from 
acetate and the repression of genes encoding enzymes unique to 
pathways of methanogenesis from methylotrophic substrates. The 
results also indicate that MreA is involved in repression of genes 
encoding the previously described Msr family of transcriptional 
regulators and MreA homologs that may also regulate genes en- 
coding enzymes unique to methylotrophic pathways. Overall, the 
results indicate that MreA plays a key role in the global regulation 
of diverse methanogenic pathways. 

RESULTS AND DISCUSSION 

Characterization of MA3302. Previous proteomics analyses re- 
vealed elevated levels in acetate-grown versus methanol-grown 
M. acetivorans (17) of the protein encoded by MA3302 that we 
designate MreA. The expression pattern of mreA was reassessed by 
two additional quantitative approaches, quantitative reverse- 
transcription PCR (qRT-PCR) and a /3-glucuronidase-based re- 
porter gene system wherein the 1,099-bp region upstream of the 
mreA translational start site (TSS) was fused to uidA. This region 
contains putative TATA and purine-rich B recognition element 
(BRE) sequences consistent with a promoter region driving tran- 



TABLE 1 Expression of mreA during growth of wild-type M. acetivorans 
strain WWM75 on acetate, methanol, or trimethylamine 



Substrate 


Expression level 




qRT-PCR (AAC r )° 


Gus (nmol-min 1 • mg protein l ) b 


Acetate 


38.0 ± 12.1 


9.70 ± 2.26 


TMA 


2.22 ± 0.68 


0.325 ± 0.075 


Methanol 


1.00 ± 0.17 


0.276 ± 0.024 



" Calculated using the threshold cycle {AAC r } method with the 16s rRNA gene used as 
an invariant control. The AC r values for methanol cells were used as the calibrator. 
Values represent the averages and standard deviations of the results of two biological 
replicate experiments assayed in triplicate. Primer and probe sequences are listed in 
Table S4 in the supplemental material. 

h /^-Glucuronidase activity was measured in cell lysate collected at mid-log growth. 
Values represent the averages and standard deviations of the results of three biological 
replicate experiments. 

scription of mreA (see Fig. SI in the supplemental material). The 
two independent approaches were in good agreement (Table 1), 
revealing expression of mreA elevated approximately 18- to 38- 
fold during growth on acetate versus methanol or trimethylamine 
(TMA). The results indicate a role for MreA specific for genes 
essential for growth with acetate. Further, the successful expres- 
sion of j8-glucuronidase validated the identification of the puta- 
tive promoter region driving transcription of mreA. 

Although MA3302 is annotated as encoding a hypothetical 
protein (24), a BLAST search revealed identity (E value, 1.9 10 -40 ) 
of MreA to the conserved domain of the COG4738 family of pre- 
dicted transcriptional regulators (25), of which none have been 
investigated. The search also revealed identity of MreA to several 
global transcriptional regulators of the archaeal TrmB (Thermo- 
coccus regulator of maltose binding) family which regulates sugar 
metabolism (2, 26) in Thermococcus kodakaraensis and Pyrococcus 
furiosus. The full-length MreA aligns with the first 111 N-terminal 
residues of the 346-residue TrmB regulator from T. kodakaraensis 
with 68% similarity, including conservation of a helix-turn-helix 
motif (see Fig. S2 in the supplemental material). Thus, MreA lacks 
the extended TrmB 235-residue C-terminal domain involved in 
sugar binding, suggesting that MreA lacks a small-molecule sensor 
domain. MreA is best described as a homolog of TrmBL2 which 
lacks the TrmB family sugar binding domain and is evolutionarily 
the most conserved regulator found in all sequenced genomes in 
the taxonomic order Thermococcales (3, 4). Although TrmBL2 of 
Pyrococcus furiosus recognizes the promoter of the maltodextrin- 
specific ABC transporter system, no pattern with respect to other 
promoter sequences can be deduced (4). Thus, it is unknown what 
DNA motifs or substrates, if any, are recognized by TrmBL2 ho- 
mologs. Remarkably, it was recently shown that the TrmBL2 ho- 
molog of T. kodakaraensis is an abundant chromosomal protein 
which binds to both coding and intergenic regions, conferring a 
thick fibrous structure contributing to the overall DNA architec- 
ture (27). Analysis of a deletion mutant showed that TrmBL2 is 
involved in global regulation, although the mechanism is un- 
known. 

Binding of MreA to promoter DNA. MreA was overexpressed 
in Escherichia coli, purified to homogeneity, and used in electro- 
phoretic mobility shift assays (Fig. 1) to determine if MreA binds 
to promoter DNA. Pronounced band shifts were observed in the 
presence of increasing amounts of MreA for DNA fragments cov- 
ering the promoter regions of mreA and for the pta and fpo oper- 
ons which represent genes from the acetotrophic and methyl- 
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Lane 1 2 3 4 5 6 

MreA 0.0 0.064 0.16 0.32 0.64 1.6 
(\iM) 




FIG 1 Electrophoretic mobility shift assays. The micromole concentrations 
of MreA are shown below the lane numbers. (A) The probe DNA encompassed 
230 bp upstream of the M. acetivorans pta start codon which included the 
TATA box, BRE, and 210 bp upstream of a putative TSS deduced from the 
previously identified TSS for Methanosarcina thermophila pta (29) which is 
conserved in M. acetivorans. (B) The probe DNA encompassed 420 bp of DNA 
upstream of the start codon for the first gene in the fpo operon which included 
355 bp upstream of the previously published TSS of the M. acetivorans fpo 
operon (36). (C) The probe DNA contained the 450-bp DNA fragment up- 
stream of the mreA start codon shown in Fig. SI in the supplemental material 
which included the TATA box, BRE, and 226 bp upstream of a putative TSS 
deduced from the previously identified TSS for Methanosarcina mazei mreA 
(35) which is conserved in M. acetivorans. (D) The probe was poly(dl-dC)/ 
poly(dl-dC). 



otrophic pathways previously shown to be regulated in response 
to the respective substrates (17, 28). Significant smearing, poten- 
tially attributable to binding of MreA at multiple sites analogous 
to abundant binding of TrmBL2 to the chromosome of T. kodak- 
araensis (27), was observed with the M. acetivorans DNA probes. 
No shift of poly(dI-dC)/poly(dI-dC) was observed (Fig.lD), con- 



firming specificity of binding to the promoter regions of M. ace- 
tivorans DNA. Addition of acetate, methanol, or TMA (10 to 
100 mM) had no significant effect on DNA binding, and manual 
inspection failed to identify a motif common to any of the pro- 
moter regions investigated, results consistent with properties re- 
ported for TrmBL2. Features of MreA which parallel TrmBL2 call 
for more-extensive investigations to determine other proteins re- 
quired to interact with MreA to sense effector molecules and de- 
termine if MreA and MreA homologs bind DNA contributing to 
chromosome architecture similarly to TrmBL2. 

Characterization of the deletion mutant strain AmreA. The 
potential role of MreA in gene expression was further investigated 
by deleting mreA in wild-type M. acetivorans strain WWM75, pro- 
ducing the AmreA strain. Growth parameters of the AmreA strain 
with methanol (Fig. 2) or TMA (not shown) were indistinguish- 
able from those of WWM75. However, a significantly reduced 
growth rate, a prolonged lag phase, and lower growth yield were 
observed for the acetate-grown AmreA strain versus acetate- 
grown WWM75 (Fig. 2), suggesting a potential regulatory role for 
MreA specific to this substrate. 

The potential regulatory function of MreA was investigated by 
comparing the genome-wide expression levels of genes between 
acetate-grown M. acetivorans wild-type strain WWM75 and the 
AmreA mutant strain by means of RNA-Seq (see Table SI in the 
supplemental material) . A total of 280 genes were identified with a 
differential of transcript abundance of >3-fold between WWM75 
and the AmreA strain (see Table S2 in the supplemental material). 
Of these 280 genes, 228 and 52 exhibited higher and lower tran- 
script abundance, respectively, in the AmreA strain and WWM75 
strains, indicating either direct or indirect roles for MreA in both 
positive and negative regulation of gene expression. A grouping of 
the 280 differentially expressed genes, several of which are essen- 
tial for either the acetotrophic or methylotrophic pathways of 
methanogenesis (see Fig. S3 in the supplemental material), into 
their respective COG functional categories revealed a dispropor- 
tionally high fraction (P < 5.85 X 10~ 12 ) in the category "energy 
production and conversion" (see Table S3 in the supplemental 
material). 

Table 2 lists genes encoding enzymes essential to the ace- 
totrophic pathway (see Fig. S3 in the supplemental material) for 
most of which the transcript abundance ratios were lower in the 
acetate-grown AmreA strain than in the WWM75 strain, indicat- 
ing that MreA is necessary for induction in M. acetivorans. Acetate 
kinase (Ack) and phosphotransacetylase (Pta) are upregulated in 
acetate-grown versus methanol-grown M. acetivorans (17), 
Methanosarcina mazei (16), and Methanosarcina thermophila 
(29), where the enzymes function in concert to activate acetate. 
The pta and ack genes are adjacent on the Methanosarcina ge- 
nomes and have been shown to be cotranscribed in M. thermo- 
phila (29). Thus, although the reduction in expression of ack was 
less than that of pta, it is likely that expression of both genes is 
dependent on MreA. Finally, the results in Table 2 indicate a role 
for MreA in regulation of carbonic anhydrase (Cam), which is 
upregulated in acetate-grown versus methanol-grown M. acetiv- 
orans (17), M. mazei (16), and M. thermophila (30, 31) and has 
been proposed to function in removal of C0 2 from the cytoplasm 
during growth on acetate. 

The reduced expression of the acetate pathway genes cam, pta, 
ack, and cdhABCDE-1 in the AmreA strain versus WWM75 alone 
(Table 2) explains the growth defect (Fig. 2) of the AmreA strain. 
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FIG 2 Growth of WWM75 and AmreA Methanosarcina acetivorans strains 
cultured with either acetate or methanol. (A) Growth with acetate. Values 
represent the averages and standard deviations of the results determined for 
five replicate cultures. (B) Growth with methanol. Values represent the aver- 
ages and standard deviations of the results determined for three replicate cul- 
tures. Symbols: 0 , wild-type strain WWM75; O, mutant strain AmreA. 



The genome of M. acetivorans is annotated with duplicate gene 
clusters (MA 1 0 1 1-MA1 0 1 6 and MA3 860 -MA3865 ) , each encod- 
ing five subunits of the Cdh complex (CdhABCDE) which is cen- 
tral to the acetate pathway. Although proteomics analysis revealed 
upregulation of both Cdh complexes (17), the complex encoded 
byMA1011-MA1016 {cdhABCDE-1) was down regulated and the 
complex encoded by MA3860-MA3865 (cdhABCDE-2) was up- 
regulated in the AmreA strain versus WWM75 (Table 2). This 
result indicates that expression of cdhABCDE-2 is independent of 
MreA and that expression is upregulated in the AmreA strain to 
compensate for decreased expression of cdhABCDE-1. Rnf and 



Mrp complexes are upregulated in acetate-grown versus 
methanol-grown M. acetivorans and proposed to be essential for 
acetate-dependent growth (28) (see Fig. S3 in the supplemental 
material); however, the results shown in Table 2 indicate no sig- 
nificant difference in expression between the AmreA strain and 
WWM75 of genes encoding either complex. It is noteworthy that 
rnf and mrp genes are absent in M. mazei and M. barkeri that 
synthesize alternative energy-converting complexes absent in 
M. acetivorans. Therefore, it is tempting to speculate that separate 
regulatory control mechanisms evolved to coordinate expression 
of the different energy-converting complexes in Methanosarcina 
species. 

There was no significant change (see Table SI in the supple- 
mental material) for the acetate-grown AmreA strain versus 
WWM75 in the expression of genes common to both the ace- 
totrophic and methylotrophic pathways (see Fig. S3 in the supple- 
mental material). This result is consistent with a role for MreA in 
the global regulation of pathway-specific genes, thereby facilitat- 
ing the switch between growth substrates. An exception was the 
gene cluster encoding the tetrahydromethanopterin 
S-methyltransferase complex (Mtr), with a 4-fold increase in tran- 
script abundance for the AmreA strain versus WWM75 (see Table 
S2 in the supplemental material) . This result was unexpected, con- 
sidering that the acetotrophic pathway has a higher demand for 
Mtr than methylotrophic pathways supported by upregulation of 
mtr in acetate-grown versus methanol-grown M. acetivorans (28). 
Thus, it appears that the role for MreA in the regulation of mtr is 
complex. 

Table 3 lists genes encoding enzymes unique to the pathways of 
methanogenesis from methylotrophic substrates (see Fig. S3 in the 
supplemental material), among which methanol-specific genes 
are upregulated in methanol-grown versus acetate grown M. ace- 
tivorans (17, 21, 22) (see Table S2 in the supplemental material) 
and M. mazei ( 16) . In the methyl transfer branch specific to meth- 
anol, the methyl group is transferred to coenzyme M catalyzed by 
MtaBC and MtaA, for which the genome contains three and two 
copies, respectively, of the encoding genes (32). In the oxidative 
branch, a portion of methyl-coenzyme M is oxidized to C0 2 , pro- 
viding electrons for reduction of the methyl group of the remain- 
ing methyl-coenzyme M to CH 4 via reactions common to all 
methanogenic pathways (33). Three genes (fmd,ftr, and mer) of 
the five encoding enzymes in the oxidative branch showed >3- 
fold increases in expression in the AmreA strain versus WWM75 
and the other two (mch and mtd) nearly 2-fold increases, indicat- 
ing involvement of MreA in repressing genes of this branch during 
growth on acetate (Table 3). Furthermore, all 14 genes encoding 
the Fpo complex involved in electron transport and energy con- 
version of all methylotrophic pathways showed increased expres- 
sion in the AmreA strain versus WWM75, supporting a role for 
MreA in repression during growth on acetate. Although proteins 
regulating the oxidative branch were not previously reported, 
methanol-specific regulators MsrA-MsrE that either activate or 
repress genes of the methyl transfer branch in M. acetivorans have 
previously been described (22). The expression of all three copies 
of mtaBC showed no significant change in the AmreA strain versus 
WWM75 strains (Table 3), indicating no direct role for MreA in 
activation of these genes, consistent with roles for MsrA-MsrE. 
No transcriptional regulators are reported for genes encoding en- 
zymes specific for methylamines in the methyl transfer branch of 
methylotrophic pathways. However, the results in Table 3 indicate 
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TABLE 2 Expression ratios of genes in the acetotrophic pathway for the acetate- 


grown AmreA strain versus the WWM75 strain of M. acetivorans 


Gene(s) 


Transcript abundance (AmreA strain/WWM75) tI 


MA0658-MA0664 (rnfXCDGEAB) 


1.13-1.33 


MA0665 (cytC) 


1.18 


MA1016-MA1011 (cdhABCDE-1) 


0.06-0.15 


MA1016 (cdhA-1) 


0.03 ± O.OO'' 


MA2536 (cam) 


0.11 (0.10 ± 0.01) b 


MA3606 (ack) 


0.4 


MA3607 (pta) 


0.18 (0.20 ± 0.06) b 


MA3860-MA3865 (cdhABCDE-2) 


2.17-4.91 


MA3860 (cdhA-2) 


5.21 ± 0.82'' 


MA4399 (cdhA-3) 


6.95 


MA4566-MA4572 (mrpABCDEFG) 


0.55-0.67 



" Determined by RNA-Seq unless noted otherwise, as indicated in footnote b. 

b Determined by qRT-PCR, calculated using the AAC r method with the 16s rRNA gene used as an invariant control. Values represent the averages and standard deviations of the 
results of two biological replicate experiments assayed in triplicate. Primer and probe sequences are listed in Table S4 in the supplemental material. Primers and probes for MA2536 
were those used previously (40). 



that MreA is involved in repression of methyltransferases specific 
for monomethylamine (MtmBC-1 and MtbA), dimethylamine 
(MtbBC-2 and MtbA), and trimethylamine (MttBC-2 and 
MtbA), for which expression of the encoding genes was greater in 
the acetate-grown AmreA strain versus WWM75. Yet genes en- 
coding duplicate methyltransferases (MtmBC-2, MtbBC-1, 
MtbBC-3, and MttBC-1) showed no significant difference in ex- 
pression levels (Table 3), indicating a different mode of regula- 
tion. Likewise, genes encoding methyltransferases (MtsD, MtsF 
and MtsH) showed no significant difference in expression levels 



for the AmreA strain versus WWM75, consistent with previously 
published roles for MsrF, MsrC, and MsrG in regulation of these 
genes in M. acetivorans (21). 

Differential expression of regulatory proteins. In addition to 
metabolic genes, several genes encoding putative regulatory pro- 
teins showed elevated expression in the acetate-grown AmreA 
strain versus WWM75, indicating that MreA is somehow involved 
in repression of these putative regulators during growth on acetate 
(see Table S2 in the supplemental material). The results suggest 
that, at least in some cases, the mreA deletion is indirectly modu- 



TABLE 3 Expression ratios of genes in methylotrophic pathways for the acetate-grown AmreA strain versus the WWM75 strain of M. acetivorans 
Gene Transcript abundance (AmreA strain/WWM75) a 

Methyl transfer branch 



MA0144-MA0145 (mtmBC-1) 


3.95-3.13 


MA0146 (mtbA) 


4.24 


MA0455-MA0456 (mtaBC-1) 


0.88-0.49 


MA0527, MA0532 (mtbBC-1) 


0.68-1.70 


MA0528-MA0529 (mttBC-1) 


0.53-0.57 


MA0859 (mtsD) 


1.36 


MA0931-MA0932 (mttBC-2) 


6.12-9.27 


MA0933-MA0934 (mtbBC-2) 


8.2-10.60 


MA1615 (mtaA2) 


0.92 


MA1616-MA1617 (mtaBC-3) 


1.21-1.22 


MA2424-MA2425 (mtbBC-3) 


0.78-0.95 


MA2971-MA2972 (mtmBC-2) 


1.37-1.44 


MA4379 (mtaA-1) 


7.83 


MA4384 (mtsF) 


1.00 


MA4558 (mtsH) 


0.93 


Oxidative branch 




MA0010 (fir) 


4.98 (7.07 ± 1. 


MA0304-MA0309 (fmdEFADCB) 


7.81-10.09 


MA0304 (fmdE) 


14.8 ± 4.0'' 


MA0309 (fmdB) 


11.3 ± 2.4'' 


MA1710 (mch) 


1.85 


MA3733 (mer) 


6.55 (8.70 ± 2. 


MA4430 (mtd) 


1.97 


Energy conversion 




MA1494-MA1507 (fpoPABCDHIJJKLMNO) 


6.35-13.64 


MA1495 (jpoA) 


9.34 ± 0.92 b 


MA1498 (fpoD) 


11.9 ± l.5 b 


MA3732 (fpoF) 


3.91 



" Determined by RNA-Seq unless noted otherwise, as indicated in footnote b. 

b Determined by qRT-PCR, calculated using the AAC T method with the 16s rRNA gene used as an invariant control. Values represent the averages and standard deviations of the 
results of two biological replicate experiments assayed in triplicate. Primer and probe sequences are listed in Table S4 in the supplemental material. 
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TABLE 4 Expression ratios of genes encoding Msr (methanol-specific 
regulatory) family and putative Mre family regulatory proteins in the 



acetate-grown AmreA strain versus the WWM75 strain of M. acetivorans 


Gene 


Transcript abundance (AmreA strain/WWM75) a 


MA0459 (msr A) 


3.06 


MA0460 (msrB) 


1.93 


MA0862 (msrF) 


8.66 


MA1671 (mreB) 


8.64 (10.2 ± 2.6) b 


MA3129 (mreC) 


8.67 (13.0 ± 3.S) b 


MA3130 (mreD) 


363.95 (1440 ± 430)'' 


MA4167 (msrH) 


11.77 (19.0 ± 5.2) b 


MA4383 (msrC) 


27.93 



" Determined by RNA-Seq (see Table S2 in the supplemental material) unless noted 
otherwise, as indicated in footnote b. 

b Determined by qRT-PCR, calculated using the AAC T method with the 16s rRNA gene 
used as an invariant control. Values represent the averages and standard deviations of 
the results of two biological replicate experiments assayed in triplicate. Primer and 
probe sequences are listed in Table S4 in the supplemental material. 



lating gene expression by affecting the expression of other regula- 
tors. 

The genes are annotated (2 1 ) as encoding nine putative sensory 
transduction histidine kinases (MA0490, MA0551, MA0619, 
MA0759, MA0863, MA1645, MA1646, MA1800, and MA2397), 
one putative response regulator (MA2445), and eight DNA- 
binding regulatory proteins (Table 4). Five of the DNA-binding 
proteins are members of the previously described methanol- 
specific regulatory (Msr) family (21, 22), consistent with a direct 
or indirect role for MreA in repression of the encoding genes dur- 
ing growth on acetate (Table 4). Locus MA0459 encodes MsrA, 
which, together with MsrB (MA0460), activates methanol-specific 
mtaBCl (22). Although msrB was elevated less than the 3.0-fold 
cutoff (Table 4), it is reasonable to expect that both expression of 
msrB and expression of msr A are directly or indirectly regulated by 
MreA since the genes are cotranscribed (22). However, the ex- 
pected increase in expression of mtaBCl by the AmreA strain was 
not observed (Table 3), indicating a requirement of other regula- 
tory factors during growth on methanol. Loci MA0862 and 
MA4383 encode MsrF and MsrC, which activate transcription of 
mtsD and mtsF, respectively, previously reported to be involved in 
the pathway of methanogenesis from dimethylsulfide (21). Thus, 
the finding that genes encoding MsrF and MsrC were upregulated 
in the AmreA strain versus WWM75 predicted increased expres- 
sion of mtsD and mtsF in the AmreA strain versus WWM75. How- 
ever, the expected increase in expression was not observed (Table 
3), indicating that the regulation of mtsD and mtsF is more com- 
plex. Expression of proteins encoded by the mtsF and mtsD genes 
is elevated 50-fold in CO-grown versus acetate-grown or 
methanol-grown M. acetivorans, suggesting a role in addition to 
metabolism of dimethylamine (36). Indeed, the protein encoded 
by mtsF (renamed cmtA) catalyzes methyl 
tetrahydromethanopterin-coenzyme M methyltransferase activ- 
ity that functions in the CO-dependent C0 2 reduction pathway 
(34). Finally, the deduced sequence of locus 4167 was found to 
have 33% to 40% identity, with conservation of the helix-turn- 
helix motif, to previously described Msr family members (21), 
identifying it as a new member of the Msr family that we designate 
MsrH. 

The three other putative DNA-binding proteins (MA 1671, 
MA3129, and MA3130) with elevated gene expression in the 



TABLE 5 Relative transcript abundance ratios for genes encoding 
putative Mre family regulatory proteins in methanol-grown or 
trimethylamine-grown versus acetate-grown wild-type M. acetivorans 
strain WWM75 



Gene 


Transcript abundance" 




MeOH'/HOAc^ 


TMA<VHOAc c 


MA1671 (mreB) 


5.09 ± 0.91 


17.7 ± 1.9 


MA3129 (mreC) 


1.89 ± 0.21 


6.51 ± 0.51 


MA3130 (mreD) 


85.7 ± 13.2 


393 ± 17 



a Determined by qRT-PCR. The AC r values for acetate-grown cells were used as the 

calibrator. Values represent the averages and standard deviations of the results of two 

biological replicate experiments assayed in triplicate. Primer and probe sequences are 

listed in Table S4 in the supplemental material. 

b Methanol grown. 

c Acetate grown. 

d Trimethylamine grown. 



AmreA strain versus WWM75 (Table 4) had deduced sequence 
identity to MreA ranging from 43% to 47% with conservation of 
the helix-turn-helix motif. Transcript abundance ratios indicated 
that each gene is upregulated in either methanol-grown or TMA- 
grown versus acetate-grown wild-type WWM75 (Table 5), a result 
in contrast to that seen with mreA, which was upregulated in 
acetate-grown versus methanol-grown or TMA-grown M. acetiv- 
orans (Table 1). The results suggest that the genes encode the 
putative homologs of MreA that we designate MreB (MA1671), 
MreC (MA3129), and MreD (MA3130) with potential to activate 
the expression of genes unique to methylotrophic pathways. De- 
spite robust increases in expression for all three homologs in the 
acetate-grown AmreA strain versus the WWM75 strain, there was 
no apparent restoration of expression of acetate-specific genes 
(Table 2), indicating that induction of the homologs in the 
acetate-grown AmreA strain does not function to compensate for 
loss of MreA. This result further supports the idea of a functional 
role for the homologs during methylotrophic growth. 

The overall results establish a role for MreA in the global reg- 
ulation of energy-converting pathways of M. acetivorans as sum- 
marized in Fig. 3. 

Conservation of MreA and homologs across the genus 
Methanosarcina. Homologs of MreA with >97% amino acid se- 
quence identity were identified in both M. mazei and M. barkeri 
(see Fig. S4 in the supplemental material). In addition, the gene 
context of mreA in both species was found to be similar to that for 
M. acetivorans strain C2A, as was the promoter DNA sequence, 
including the putative TATA box and BRE previously mapped for 
M. mazei mreA (35) (see Fig. SI in the supplemental material). 
Together, these data are consistent with similar functions for 
mreA in all three species. M. barkeri has a second MreA homolog 
designated MreA-2 with 88% sequence identity to MreA of 
M. acetivorans (see Fig. S4 in the supplemental material). How- 
ever, unlike mreA-1, the promoter region and gene context for 
M. barkeri mreA-2 were found to be dissimilar to those of the other 
mreA homologs. Thus, it is unclear whether the functions for 
M. barkeri mreA-1 and mreA-2 are similar. MreB, MreC, and 
MreD are also conserved in both M. mazei and M. barkeri with 
>80% identity. Conservation of MreA, MreB, MreC, and MreD 
indicates an important role for the Mre family consistent with 
regulation of energy metabolism across the genus Methanosar- 
cina. 
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FIG 3 Global regulatory network model of energy-converting pathways. Genes 
mreB, mreC, and tnreD (Methanosarcina regulator of energy-converting metabo- 
lism) are homologs of tnreA postulated to encode proteins involved in regulation 
of genes required for conversion of methylamfnes and methanol to methane. En- 
zymes for conversion of acetate to methane are encoded by ack (acetate kinase), pta 
(phosphotransacetylase), cdh-1 (carbon monoxide dehydrogenase), cam (car- 
bonic anhydrase), and tntr (membrane-bound tetrahydromethanopterin 
S-methyltransferase complex) . Enzymes required for conversion of methylamines 
to methane are encoded by mtbA (methylamine-specific methylcobamide:CoM 
methyltransferase), mtmBC-1 (B, monomethylamine-specific methyltransferase; 
C, cognate corrinoid protein), mtbBC-2 (B, dimethylamine-specific methyltrans- 
ferase; C, cognate corrinoid protein), and mttBC-2 (B, trimethylamine-specific 
methyltransferase; C, cognate corrinoid protein). Enzymes required for conver- 
sion of methanol to methane are encoded by mtaBC-1, mtaBC-2, mtaBC-3 (B, 
methanol-specific methyltransferase; C, cognate corrinoid protein), and mtaA-1 
(methanol-specific methylcobamideiCoM methyltransferase). Enzymes common 
to the pathways for conversion of methanol or methylamines to methane are 
encoded by mer (F^o-dependent methylene-H 4 MPT reductase), mtd (F 420 - 
dependent methylene-H 4 MPT dehydrogenase), mch (mefhenyl-H 4 MPT cyclohy- 
drolase), fir (formyhnethanofuran; H 4 MPT formyltransferase), fmd 
(molybdenum-dependent formyl-methanofuran dehydrogenase), fpo (coenzyme 
F 420 dehydrogenase complex), zndfpoF (a subunit of the coenzyme F 420 dehydro- 
genase complex). Proteins regulating genes specific to the conversion of methanol 
are encoded by msrA, msrB, msrD, msrE, msrC, and msrF (methanol-specific reg- 
ulators). An enzyme proposed to be specific for conversion of dimefhylsulfide to 
methane is encoded by mtsD (methyltransferase). An enzyme in the pathway for 
conversion of carbon monoxide to methane is encoded by cmtA (cytoplasmic 
methyltransferase) (formerly mtsF). 



Conclusions. The results presented here have identified a 
novel archaeal regulatory protein (MreA) which was previously 
annotated as a hypothetical protein. MreA is the first global regu- 
latory protein described for methanogenic species. Overall, the 
results presented are consistent with MreA functioning in the 
global regulation of energy-converting metabolism based on the 
finding that MreA is involved in regulation of (i) genes unique to 
the acetotrophic and methylotrophic pathways, (ii) Msr family 
regulatory proteins, and (iii) homologs of MreA that are upregu- 
lated in methanol-grown or TMA-grown versus acetate-grown 
M. acetivorans (Fig. 3). 

Comparative transcriptomic analyses revealed a dual role for 
MreA in both the transcriptional activation of genes specific for 
the metabolism of acetate and repression of genes unique to 
methylotrophic pathways. Furthermore, the results indicate 
MreA is involved in the transcriptional repression of genes encod- 
ing Msr family regulators and MreA homologs with the potential 
for regulating genes unique to methylotrophic pathways. How- 
ever, a direct versus indirect role for MreA in regulating gene 
expression has yet to be determined. MreA joins a growing num- 
ber of archaeal proteins, including TrmBLl, Tgr, and SurR (2, 3, 
5), that function in both the repression and activation of gene 
expression for which the mechanism has not been reported. The 
similarity of MreA to TrmBL2 warrants more-extensive investiga- 
tions to determine other proteins required to interact with MreA 
to sense effector molecules and the manner in which MreA and 
homologs bind DNA contributing to chromosome architecture 
similar to that of TrmBL2. 

MATERIALS AND METHODS 

Cell growth. Growth and harvesting of M. acetivorans C2A (DSM 800) in 
a single-cell morphology cultured with high-salt (HS) broth medium with 
acetate, methanol, or TMA were performed as previously described (36). 
Growth of all strains was monitored by measuring the optical density at 
600 nm. 

RNA isolation. All M. acetivorans RNA samples were isolated with an 
RNeasy Total RNA Minikit (Qiagen). Purified RNA was treated twice with 
RNase-free DNase I (Qiagen) and once with RQ1 DNase (Promega) to 
remove contaminating DNA. 

TaqMan quantitative RT-PCR. TaqMan assays were performed as 
previously described (36) with total RNA isolated from the mid- 
exponential phase (A 600 ~ 0.15 to 0.2, 0.4 to 0.5, or 0.4 to 0.5 for acetate-, 
TMA-, or methanol-grown cells, respectively). Primer and probe se- 
quences are listed in Table S4 in the supplemental material. 

Construction of the M. acetivorans mutant strain. Liposome- 
mediated transformation and homologous recombination-mediated 
gene replacement were performed as previously described (37, 38) to gen- 
erate an M. acetivorans AmreA::pac-hpt strain (here designated the AmreA 
strain) and VmreA-uidA. The AmreA mutant strain was generated in an 
M. acetivorans WWM75 background, and a VmreA-uidA promoter fusion 
strain was constructed in an M. acetivorans WWM73 background. Cells 
were transformed with 2 /ng appropriate shuttle vector DNA (see Table S4 
in the supplemental material) (linearized by restriction digestion in the 
case of the gene knockout shuttle vector) , and transformants were selected 
on HS agar media with 1.5% (wt/vol) agar, with 50 mM TMA and 
2 jug ml -1 puromycin (Sigma, St. Louis, MO) added from sterile, anaer- 
obic 100X stock solution. 8-Aza-2,6-diaminopurine (8ADP) (Sigma, St. 
Louis, MO) sensitivity was tested in media containing 20 jug ml -1 8ADP 
added from sterile, anaerobic 200 X stock solution. Removal of pac-hpt 
cassettes in mutant strains by markerless exchange was not performed in 
this study. 

Gus assays. /3-glucuronidase (Gus) assays were performed essentially 
as previously described (23, 39, 40). An M. acetivorans strain carrying the 
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PmreA-uidA reporter gene was grown to mid-exponential phase on dif- 
ferent methanogenic substrates (A 600 ~ 0.15 to 0.2, 0.4 to 0.5, or 0.4 to 0.5 
for acetate-, TMA-, or methanol-grown cells, respectively) and harvested 
by centrifugation for 10 min at 5,000 X g. Cells were osmotically lysed by 
addition of 50 mM Tris-HCl (pH 8.0) containing 1 mM dithiothreitol 
(DTT) and 0. 1 fxg mH DNase. Cell debris was cleared by recentrifugation 
for 15 min at 16,000 X g. Cleared lysate was held at 4°C or frozen at — 80°C 
until use. For each assay, cleared lysate was diluted 1:20 into 50 mM 
Tris-Cl (pH 8.0) prechilled to 4°C and brought to room temperature. 
Assays were initiated by addition of 4 mM p-nitrophenyl-j8-d- 
glucuronide, and activity was monitored by measuring absorbance at 
415 nm with a Beckman DU 640 spectrophotometer. The specific activity 
(in nanomoles per minute per milligram of protein) was based on a pre- 
viously reported molar extinction coefficient (39). The protein concen- 
tration of the cell lysate was determined using a bicinchoninic acid (BCA) 
protein assay kit (Pierce) with bovine serum albumin (BSA) as the stan- 
dard. 

RNA-Seq. Total RNA from acetate-grown wild-type and \mreA 
M. acetivorans was collected at the mid-exponential phase (A 600 ~ 0.15 to 
0.2). Transcriptome libraries were prepared by following the Whole Tran- 
scriptome Library Preparation for SOLiD Sequencing Protocol (Rev. B) (Ap- 
plied Biosystems). The libraries were constructed with reagents from a 
SOLiD small RNA expression kit (Applied Biosystems) . Briefly, total RNA 
samples (1 ju,g) were fragmented with RNase and purified. Fragmented 
RNA (100 ng) was ligated overnight to adaptors and subsequently reverse 
transcribed. Purified cDNA was size selected and then amplified from 
bar-coded primers. The resulting material was purified and assessed for 
proper amplification using an Agilent Bioanalyzer DNA 1000 chip (Agi- 
lent Technologies). Bar-coded libraries were combined and prepared for 
sequencing according to the Applied Biosystems SOLiD 3 Plus System P em- 
plated Bead Preparation Guide (Applied Biosystems). Fragment library 
sequencing (50 bp) was performed on a SOLiD 3 Plus system according to 
the Applied Biosystems SOLiD 3 Plus System Instrument Operation Guide 
(Applied Biosystems). Sequenced library fragments were mapped to the 
M. acetivorans reference genome sequence with the Qseq program in Ar- 
rayStar 4.0. Noncoding RNAs (e.g., rRNA, tRNA, miscellaneous RNAs) as 
well as nonunique reads were excluded from mapping. Expression values 
for each gene were normalized by RPKM (assigned reads per kilobase of 
feature per million mapped reads). Statistical analyses of the normalized 
data were performed using a moderated t test with false-discovery-rate 
(FDR) multiple-testing correction (Benjamini-Hochberg) to determine 
differential transcript abundance. Changes in transcript abundance were 
considered significant if they were >3-fold (P values < 0.05). Functional 
categories represented by genes exhibiting differential transcript abun- 
dances in the AmreA strain were evaluated using a hypergeometric distri- 
bution to determine if any particular categories were statistically over- or 
underrepresented. 

There were 4,542,508 and 2,099,564 unique reads assigned for the two 
respective biological replicates of the wild-type strain and 2,704,850 and 
2,941,777 unique reads assigned for the two respective biological repli- 
cates for the mutant strain. 

Cloning, expression, and purification of MreA. The gene encoding 
MreA was amplified from M. acetivorans genomic DNA by PCR. The 
PCR-amplified DNA fragment was cloned into pSUMO vector (LifeSen- 
sors, Malvern, PA), generating plasmid pSUMO-MreA. pSUMO-MreA 
contains the N-terminally His-tagged Sumo domain-MreA fusion. 

The His 6 -Sumo-MreA fusion was overproduced in E. coli BL21(DE3) 
cells transformed with pSUMO-MreA. Cells were grown in LB at 37°C 
with shaking at 250 rpm until an optical density at 600 nm of 0.5 was 
reached, after which the culture was induced with 1 mM isopropyl-j3-d- 
thiogalactopyranoside (IPTG) and then harvested by centrifugation 16 h 
after induction. Approximately 5 g (wet weight) of cells was suspended in 
15 ml of 20 mM Tris-Cl (pH 8) containing 150 mM NaCl, 20 mM imida- 
zole, 1 mM phenylmethylsulfonyl fluoride (PMSF), 20 u,g ml -1 DNase, 
and 10% (wt/vol) glycerol (buffer A). The cells were lysed by two passages 



through a French pressure cell at 138 MPa. The lysate was centrifuged at 
74,000 X g for 30 min at 4°C. The supernatant solution containing the 
His 6 -Sumo-MreA fusion protein was filtered (pore size, 0.45 /j,m) and 
applied at a flow rate of 0.5 ml/min to a column containing 25 ml of Ni 
Sepharose 6 fast-flow resin (GE Healthcare). The column was then 
washed with 200 ml of buffer A a flow rate of 2 ml/min. His 6 -Sumo-MreA 
was then eluted from the column with 60 ml of 20 mM Tris-Cl (pH 8) 
containing 150 mM NaCl, 500 mM imidazole, and 10% (wt/vol) glycerol 
(buffer B). The eluate was concentrated to 2.5 ml with a Vivacell concen- 
trator fitted with a 10,000-molecular-weight-cutoff filter. The concen- 
trated His 6 -Sumo-MreA was desalted with a PD-10 gel filtration column 
(Amersham Biosciences) developed with buffer A. MreA was cleaved 
from His 6 -Sumo by incubating the full-length fusion protein at 30°C for 
2 h with Sumo protease 1 (LifeSensors, Malvern, PA) (10 units added per 
100 u,g protein) plus 2 mM DTT. The cleaved protein was applied at a flow 
rate of 0.5 ml/min to a second column containing 25 ml of Ni Sepharose 6 
fast-flow resin (New England Biolabs) . The flowthrough containing MreA 
was collected and concentrated to 2.5 ml with a Vivacell concentrator 
fitted with a 10,000-molecular-weight-cutoff filter. The potein concentra- 
tion was determined with a BCA protein assay kit (Pierce) and bovine 
serum albumin as the standard. 

Electrophoretic mobility shift assay. Probe DNA (100 nM final con- 
centration) was titrated with the indicated increasing concentrations of 
MreA in 10 u.1 binding buffer (10 mM Tris-Cl [pH 8.0], 100 mM NaCl, 
0.1 mM EDTA, 5 mM DTT, 15% glycerol) containing 0.6 jug bovine 
serum albumin (BSA) and 2.5 u,g heparin. Poly(dI-dC)/poly(dI-dC) 
(Sigma-Aldrich) was used at 100 nM (final concentration). Reaction mix- 
tures were incubated 30 min at 21°C and loaded onto prerun 6% native 
polyacrylamide gel electrophoresis (PAGE) gels containing 15% glycerol. 
Gels were developed with 0.5 X TBE buffer (45 mM Tris-Cl, 45 mM bo- 
rate, 1 mM EDTA, 15% glycerol) until the desired migration was achieved. 
Following electrophoresis, gels were stained for 20 min in 1 X SYBR green 
stain diluted in 0.5 X TBE. Gels were visualized by excitation at 254 nm 
with a UV epi-illumination lamp and a 490-nm Longpass SYBR photo- 
graphic filter. 
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